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Abstract 

Global existence to the coupled Einstein-Maxwell system which rules the dynamics of a kind 
of charged matter with a pseudo-tensor of pressure is proved, in Bianchi I-VIII spacetimes. We 
study the geodesics completeness, the asymptotic behavior, the positivity conditions, and we 
prove that the problem is well-posed in the sense of Hadamard. 
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1 Introduction 

In relativistic kinetic theory, global dynamics of several kinds of charged and uncharged matter 
remain an active research domain in General Relativity, in which cosmological constant plays a 
central role for astrophysical reasons. In the present paper, we consider the Einstein equations 
with in the sources, the pseudo-tensor of pressure due to A.LIGHNEROWIGZ [5] which is the 
general form of a relativistic fluid tensor. We prove the local and global existence and we study 
the asymptotic behavior. We study the positivity conditions and we prove that the problem 
is well-posed in the sense of Hadamard, which means that the global solution is a continuous 
function of the initial data. The work is organized as follows: 

In section 2, we introduce the problem and we give the evolution and the constraints systems. 
In section 3, we study the constraints, dehne the initial values problem and introduce the 
relative norms. In section 4, we construct the local solution by iterated method. In section 5, 
we prove the global existence theorem. In section 6, we study the asymptotic behavior and the 
geodesic completeness. In section 7, we study the positivity conditions. In section 8, we study 
the well-posedness in the sense of Hadamard. 
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2 The problem.Evolution and constraints systems 


2.1 The problem 

We are on Bianchi, time oriented spacetimes, which are (M, '^g) spacetimes of type I-VIII, where 
is the metric with signature (—, -|-, -|-, -|-) and M = M x G where G is a three dimensional 
connected Lie group, "^g has the form; 

"^g = —dt^ + gij{t)e^ ® ( 2 . 1 ) 

where (cj) is a left invariant basis in G and (e*) the dual basis; g = {gij) is a positive dehnite 
3-dimensional metric on G and we adopt the Einstein summation convention A°‘Ba = E A°‘Ba- 

a 

The Greek indexes a, 13, ... range from 0 to 3 and the Latin indexes i, j, ... from 1 to 3. The 
index 4 {as "^g) is for quantities on M. The vector n = dt being orthogonal to G, we complete 
the basis (e^) on G to obtain a basis (n, e^) = (clo, ^i) on M with; 

^eo = dt ; = e^; = *e^dx (2.2) 


where 

= 1; = 0; = e} (2.3) 

The structure constants of the Lie algebra Q of the Lie group G are denoted and dehned 
by: 


[G,e,] = Gj.e, (2.4) 

where [ , ] is the Lie brackets in Q. Due to the antisymmetry of [ , ], we have: 

Ci = -Gj) (2.5) 

If is the covariant derivative in ^g, the Ricci rotation coefficients dehned by: 

( 2 . 6 ) 


If and ^Tq, are tensors on M, we have: 

= ^e„( 

= 4e„( ^Tx 

We study a huid, subject to the following system in which t is the only independent variable; 


(2.7) 
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‘^Rap ~\'^R ‘^9ap + A ^Qap = 8TT{‘^Tap + (2.8) 

^ = e V, (2.9) 

^Vp + "V., = 0; (2.10) 

where; 

• (2.8) are the Einstein equations in ( *gap) with ^Rap the Ricci tensor, ‘^Rap 

the scalar curvature, A a constant called the cosmological constant, ^Tap the Maxwell tensor 
associated to the electromagnetic held ^F = ( ^Fap) and dehned by: 

"Fa^ + ^F„a "F/ (2.11) 

and ^Tap the symmetric 2-tensor; 

‘^Tc.p = + '^Qcp ( 2 . 12 ) 

where p > 0 is an unknown function of f, called the proper density, ^-u = ( ^u") is the unknown 
material velocity of the huid, is a unit vector oriented towards the future direction, and ^Qap 
is a symmetric 2-tensor called the pseudo — tensor of pressure, due to A.LICHNEROWICZ 
[5]. ^0Q,p = 0 is the case corresponding to pure matter, and ^0Q,p = p ^Qap where p is a scalar 
function representing the pressure, corresponds to a perfect relativistic huid. We adopt on ^0Q,p 
the assumptions: 

"^0"^ = p ^u^ and '‘0oa = ‘^Zq ^ (2-13) 

where ( ‘^Za) is a future pointing vector, <5° is the Kroneker’s symbol and is a similar symbol 
such that = 1 and = 0. The hrst formula in (2.13) is due to A.LICHNEROWICZ 
and the second formula which gives "^©oo = ( ^Zq^ and hence, "^Too = |p ( ^Ug)^ -|- ^©gg > 0 
will be very helpful. We suppose that: 

^9ij (2.14) 

from where, we have: 

4 1 

4rp ^44 I 4 

Tij = -p Ui Uj + -p Pij 

which is the stress-tensor of a perfect huid of pure radiation. 

• (2.9) and (2.10) are the hrst and second group of the Maxwell equations, for the electro¬ 
magnetic held ^F = ( ^Fjj) which is an antisymmetric closed 2-form. ^F°* and ^Fjj are 
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respectively, the electric and magnetic parts of 

In equations (2.9), e > 0 is an unknown function called the Maxwell current created by the 
charged particles. 

(2.10) expresses only the fact that d = 0. 

Let us recall that to solve the system (2.8)-(2.9)-(2.10) is to determine all the unknown functions 
^Fai 3 , Pi "^0a/3 and 6, which depend only on the time variable t. 


2.2 The equations 

The unit vector satishes the relation; 


V = -1 (2.15) 

which gives, since is future oriented and using (2.1); 

V = v^l + Qij (2.16) 

We study the Einstein equations in 3+1 formnlation, which means that, they are seen as giving 
the evolution of the triplet (Xlti do-^t); where S* = {f} x G, gt = {gij(t)) is called the hrst 
fundamental form of S*, and Kt = {Kij{t)) is called the second fundamental form of In the 
present case, Kij is dehned by; 

= -\9tgij- (2.17) 

We now introduce a very useful quantity called the mean curvature and dehned by; 


H = 


(2,18) 


Let us give the expressions of as they are in [1]. They are 


^7oo= 

4 / 1 

hj ■ lij 


We deduce from (2.19) that; 


= 0 ; ^ 7 ° = -K^r, Ho = H 

7“ [-+,*». + . 

+ . = CL- - V-- = C- 

111 hi 111 ^11' 


= -K^- 


(2.19) 


( 2 , 20 ) 
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2.2.1 Equations in p, and e 

We always have: 

Hence (2.8) implies the conservation conditions; 

= 0 


( 2 . 21 ) 


from where we deduce, using expression (2.12) of and the assumptions (2.13) on ^0"^: 


+ - ^V«(p + = 0. 

o 


( 2 . 22 ) 


Now a direct calculation gives: 

(2.23) 

^Va{p ^Va{p + (p ^V«( V). (2.24) 

We deduce from (2.21), using, (2.23)-(2.24) and the Maxwell equation (2.9): 


e + q [ ^V„(p "‘u") -h (p ^u") '‘Vq,( -h p = 0. 


(2.25) 


The contracted mnltiplication of (2.25) by gives, using ‘^F^y^ = "^Fpx = 0 

(since is antisymmetric) and ^-up = —1 which implies ^up ^Vq, = 0; 


V«(p 4u“) + p = 0 


If we return to (2.25), we obtain; 


I (p + e ^F\ = 0 


We deduce from (2.26) that; 


(2.26) 


(2.27) 


p = A{t)p 

with A{t) = - ( ^ ^ 




-k: + q 




4^0 


4 4^0 4.^0 

where the dot denotes the derivative with respect to t. Integrating (2.28), we obtain; 


p = p(0) exp ( I A{s)ds 


(2.28) 

(2.29) 
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which shows that: (p = 0) (p(0) = 0). In what follows, we suppose that; 


p(0) > 0 


(2.30) 


which implies p > 0. (2.27) gives: 


V V„( V) = (2.31) 

which is the differential system of current lines. It shows that in the vacuum ( = 0), 

satishes the geodesics equation; ^Vq,( = 0. Equation (2.31) gives, taking (3 = 0 and 

(3 = using (2.19) and (2.7); 


4.-,0 


= K 




W 364 qV 
^ 4p ' ^ 

, 3 6 


F = 2K] 




+ 


4p 




3e 

4p 


r-! I ^ 

j 4^0 


(2.32) 

(2.33) 


and equation (2.26) in p writes; 








-Kl + Cl,j^]p--e^F, 


3 4^0 V 


(4^0 


,2 • 


(2.34) 


Next, by (2.14) we have: ^0*-^ = | so, is given by p and We then only look for 
40O«, We deduce from (2.13), (2.7), that "^0°" satisfy the system; 


I 4000 _ j^i 4000 _ 4 q 0 j ^ \pgiJKij + p 4^0. ^2.35) 

I 400* ^ ^ 400* ^ 2ir* "0°^ - + 7;,p^'') + p ^U\ (2.36) 

Now the electromagnetic held always satishes the identities ^F"^ = 0. This implies, 

using the Maxwell equation (2.9): 

(e S“) = 0. (2.37) 


(2.37) gives: 

e = - (i^ "" 

which integrates to give: 


e{t) = e(0) exp 




4^0 



This shows that: 

(e(0) > 0) 77 (e > 0). 


(2.38) 


(2.39) 
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In what follows we adopt; 

Co := e( 0 ) > 0 . 

Using the equations (2.32)-(2.33) in ^^ 1 “, equation (2.38) in e can write; 

3 


e = 


A F 

u 


yPY + ( 4 ^ 0 ) 2 ^ • 


(2.40) 


(2.41) 


2.2.2 Equations in ^Fij and constraints 

The Maxwell equations (2.9) for /3 = i, namely ^Vq, = e can write using (2.7) and 
(2.19), to give the equation for the electric part; 

4^0i ^ ApOi _ Apki _ 1 ipjk _ 

Now, the Maxwell equations (2.10) split into the equations; 

n 

JO 


Vo + V, ^Fp + "V, Vo, = 0; 




(2.42) 


V, ^Fp + V, V,, + Vfc V,, = 0. 

(2.43) can write using (2.7), and (2.19), to give the equation for the magnetic part; 

% = V°'. 

Next, for /9 = 0, the Maxwell equations "^Vq, = e * vP gives the constraints; 


(2.43) 

(2.44) 

(2.45) 


Ci, + e = 0 (2.46) 

and (2.44) gives the constraints; 

CP ^F,i + *Fa + Cl, V,, ^ Cf,^ = 0. (2.47) 

2.2.3 The Einstein equations in 3+1 formulation. Notations.System 

The principle of the 3+1 formulation is to deduce from the Einstein equations which are a 
second order partial differential equations system, an equivalent first order system in gij and 
Kij. In the uncharged case, f^F = 0), the calculation is classical. We did the charged case 
(^F 7 ^ 0) in [ 6 ] and it is adopted here without difficulty. So are the constraints. Expressing the 
tensor we obtain; 

^ Voo = VO* vo^' + 1(7*" +' Vfc; V,, 

< Voi = - VO" V,fc (2.48) 

4^ {i T^Ok 4 pOZ 1 „ „km„nl 4 rp 4 jp i „kl 4 jp 4 jp 

Fj — CidijQkl — Qikgjl) r r — jdij 9 9 ^ki rran+9 Fik rji 
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For we have: 

4 4 4 1 

^00 3^ ^00? T{)j uq Uj Tij Ui uj ^p pij. 


(2.49) 


Since the indexes are now clearly specihed, we will denote in what follows; 


4^0* 

= Fb 

^F - = 

j- ij 

= Fij] ^Tqo — 

. An^ _ . Arp _ rp 

00) -^02 02) ^ ij ij 

^Do 

= T'oo; 

4 

To* = 


Rji 

V: 

= vP] 

= 

F 


4000 

= 000 

; ^00* 

= 00*; 400 

= 0*-?' 

*9i3 -- 

= gij 





and we will have for and F)®, using (2.50); 


(2.50) 


I = g^^g^^Fim 

\ 4f,o = = ^g,, 

Now using expression (2.48) of (2.49) of the notations (2.50), equations (2.33) of 
equation (2.34) of p, equations (2.35)-(2.36) of equation (2.41) of e, equation (2.42) 
and (2.45) of ^F°* and ^Fjj, the constraints (2.46) and (2.47), reference [6], and the notations 
(2.50), we obtain the evolution system (S); 


{S) 


kij = R, 


'U 


g,, = -2F,,; (2.51) 

HKij — 2KjKii — 87r{Tij + Tij) + 47r(—F qo + g^^Tim)gij — Agij',{2.52) 

(2.53) 

(2.54) 


4 - = Cf^guk; 


• 7/^ • 1 

E‘ = HE' - Cj,E‘’- - - -qtg’'g'‘'"F,„; 


F = 2KjU^ — 7 


u 


3 1 




a^E^g^^F^iU 


uo 4pu0 jk 4p(u0^2 jk 


,,FR. 

5 


000 _ ^000 _ Q 0Oi ^ ^ ^^0. 

00* ^ ^00* ^ 2Fi0O^' - + pF; 

3 ip ')p ip 3 iP 


e = — 


KiiEu^ , 


012 


■ + 




, 3 FE^ 

-K,]e + 


(2.55) 

(2.56) 

(2.57) 

(2.58) 

(2.59) 


Where in (2.52), Rij is the Ricci tensor associated to gij and whose expression due to R.T. 
JANTZEN [3] is: 

(2.60) 


kj = iLr 


pi ijL imi 


-7,™. 

mj Hi 
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Next, we obtain the constraints; 

Ft — KijK^^ -\- F[^ = 167r(Too + ^oo) “1“ 2A, (2.61) 

= -87r(ro,+To,), (2.62) 

< 

CljFki + Cji^Fii + Cl^Fji = CjjjFfcp = 0, (2.63) 

, Cl,F^^ + eu° = 0, (2.64) 

where R = g^^Rij and V is the Levi-Civita connection associated to 5 ' = {gij)- The constraint 
(2.61) is called the Hamiltonian constraint. As we will see, this constraint is fundamental. 

3 Study of constraints.The initial values problem.Relative 

Norms 

3.1 Study of constraints 

We prove that, if the constraints are satished by the solutions of the evolution system at t = 0, 
then the constraints are satished everywhere. Let us set: 


A = R — KijK^^ -\- — 167r(Too + log) — 2A, 

Aj = W^Kij + 87r(roj -|- T^j), 

■* Aijk = 

B = Cl^F^^ + eu^, 

^ W = (A, A,, A,,, By. 

Proposition 3.1. 1) The quantities A ,Aj: Aijk, B satisfy the relations: 


' A = 2HA+ 2g^A-jA,, 

( 3 , 1 ) 

Aj = F[Aj, 

( 3 , 2 ) 

Aijk 0 , 

( 3 , 3 ) 

, B = HB. 

( 3 , 4 ) 


2) The solutions of the evolution system satisfy the constraints everywhere, if and only if they 
satisfy the constraints at t = 0. 
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Proof. 1) See Appendix . 

2) If we set W = (A, Aj, Aijk, By, then W satishes: W = LW, where L is a matrix. So we 
have: 

W{t) = W{0) exp L{s)ds 

This shows that {W = 0) -v^ (ID(0) = 0). This means that the constraints are satished 
everywhere if and only if they are satished at t = 0. In what follows, we consider that the 
constraints are properties of the solutions of the evolution system. 

□ 



3.2 The initial values problem 

Let; = (AT® ), = (F°) be given 3x3 matrices with g^ positive dehnite, 

symmetric and antisymmetric matrices. 

Let: 0^ = (6'°’") be given vectors. 

Let: /9° > 0, and e° > 0 be given numbers. 

We look for solutions; g = {gij) positive dehnite 3x3 matrix, K = {Kij) symmetric 3x3 
matrix, F = {Fij) antisymmetric 3x3 matrix, u = (u*), E = (i?*), 0 = (0°“) vectors, p > 0, 
and e > 0 two functions, such that at f = 0, we have: 

I p( 0 ) = gy K( 0 ) = Ky F( 0 ) = Fy Eiy) = Ey 

y m(0 ) = [/O; 0(0) = 0°; p(0) = pO; e(0) = 6°; 

For t G [0,T[, T < -|-cxd. We will set: 

[/o-o = y/f + c/O [/0-* f/O-T = v^l + gijU^u^ (3.5) 

K^, F^, E^, U^, 9^, p° and e° are the initial data. In what follows, we suppose that they 
satisfy the constraints at f = 0. There are eight initial data for the four constraints. This means 
that there are four degrees of liberty to the initial data. 

3.3 Relative norms 

We now introduce the notion of relative norms due to RENDALL [7]. 
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Lemma 3.1. Define the norm of the n x n matrix A by: 

||yl|| = ^ 0,a; e . 

If Ai and A 2 are two nxn matrices, Ai positive definite, the norm of A 2 with respect to Ai is: 

||dl 2 |Ui = ^ 0 ,a; e . 

From the dehnition we have: 


P2||<P2|U,II-4 


( 3 ^ 6 ) 


We also have: 

P 2 IU 1 < [Tr{Af^A2Af^A2)Y (3.7) 

If yl = {oij) is a n X m matrix, one dehnes another norm by: 


lAI = sup{\aij\-,i = 1,2, = 1,2, . (3.8) 

Lemma 3.2. Let (m“) = where vfi and u* are linked by (3.5). Let be given. Then, 

there exists a constant C > 0 such that: 

7 /^ 3 1 3 

- <^1^1^; \F^^\<{grsF^^F^^)-^\g\-^. (3.9) 

u'^ I I \ / 

Proof. Take in lemma 3.1, Ai = A 2 = (a*-^) with a** = u* and = 0 if f 7 ^ j. A direct 
calculation, using Af^ = {gifi gives: 

Tr{Af^A2Af^A2) = (3.10) 


(3.6) and (3.7) then give, using (3.10): 

II2I2II < (a“o«)’ llsl 


( 3 . 11 ) 


But we have: 


(a*%ij)5 = (a**an)5 

= (fl-la^a**)^ < \g\^ [gikFu^Y . 

The hrst relation (3.9) is due to (3.11) using the inequalities |u*| < IA 2 I < (711^211, since all 
the norms on a hnite dimensional vector space are equivalent, the relation = (l -|- gi^Fu^) ^ 
and the fact that \g\ and || 5 r|| are equivalent. For the second relation (3.9), take a** = and 
a*-^ = 0 if i ^ j and proceed as above. This ends the proof of lemma 3.2. □ 
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4 Local existence by iterated method 


4.1 Construction of the iterated sequence 


Consider the initial values in paragraph 3.2. Set: 


I goit) = g^- Koit) = Fo{t) = Eo{t) = E^- 

y Uoit) = 0o(t) = 0^] poit) = po; eo(t) = cq; 

We also set: Vq = {go, Kq, Fq, Eo,uo,Qo, Po,eo). For n e N, ifl4 = ©n, Pn, e 

is known, then dehne W+i by: 

K+l = Vo+ f f{Vn{s))ds 
Jo 

where f{Vn) is the right hand side of the evolution system (S) in which: 

*) < = + gn,ijKui > 1 ; 

*) 7n,ii is obtained by replacing in (2.19) gij by gn^ij] 

*) Rn,ij, Tn,a/ 3 , Qn,ij are dehned by the same method. 

We then obtain the iterated sequence (W) which is of class on a maximal interval [0,T„[, 
> 0 where = {gn,ij) is symmetric and positive dehnite, iC„ = {Kn^ij) is symmetric, 
Fn = {Fn,ij) is antisymmetric, Pn > 0 and > 0. 

4.2 Estimation of the iterated sequence 

Proposition 4.1. There exists a number T > 0, independent of n, such that, the iterated 
sequence W = (Pn, Kn, Fn, En, Un, ©n, Pn, Gn) is defined and uniformly bounded on [0, T[. 

Proof. Let N eN*. Suppose that for n < iV — 1, we have the following inequalities: 

^ |p„-p°| < di; \Kn-K^\ < d2; |Fn - F^l < dg; 

< \En-E^\ < d4; \un-U^\ < ds; |0n-0°| < de; (4.1) 

|en-eo| < Af, |pn-po| < dg; {detgn)~^ < dg 
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where Ai,i = 1, 2, 3,4, 5, 6, 7, 8, 9, are strictly positive constants. We prove that we can 
choose these constants such that we still have the inequalities (4.1) for n = N, for T sufficiently 

small. 

Note that for t sufficiently small, we have 

|pv-i-p°|<4 

from where we deduce that 

-A— < A 

PN-l — P° ■ 

we take Ag = Taking the inequalities (4.1) into account, the definition of the iterated 
sequence, the expressions of 

Hn-i-, RN-i,ij, T]\f_i^ij, and u%_i = \J^ 

we can find constants Sj, i = 1,2, 3,4, 5, 6, 7, 8 strictly positive depending only on Ai such that; 


1 \9N,ij\ < 

Bi] 

1 1 

< B 2 ] 

1 

1 < 

B 3 ; \E}, 

< Ba] 

(4.2) 

VI 

B,] 

|e»l 

< Be] 

ev 

< 

By] \pn\ 

< Be- 


By integration of (4.2) 

we have: 







\9n - 9^\ < Bit] 

\Kn 

-K^\ 

< B2t] 

\Fn- 

F0| 

< Bet] 

\En-E^\ 

< BA] 

\un-U^\ < B^t] 


-e^\ 

< Bet] 

Cv — Co 

< BA] 

\Pn - Pol 

< BA- 


(4.3) 


Let us bound (det^^Ar) By the definition of the iterated sequence we have: 

d 

= —2^KM-i,ij- (4.4) 

On one hand we have the equality: 

^[In(det^jv)] = 9%j^9N,ij- 

On the other hand we have: 

^[ln(det5(jv)] = (det ^Ar)"^^(det ^fTv) 

= -det5fAr^(det^Ar)"L 

The relation (4.4) then implies: 

^(detc/Af)"^ = {2g'^KN-i,ij)idet gN)~^ 
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which is a differential equation in (det g^) ^ on [0; t], t > 0, whose solution is; 

(det^fTv)"^ = (det5f°)"^exp 2g"^KN-i,ij{s)ds^ . (4.5) 

The relation (4.4) which is similar to (2.17), shows that g^ and are the hrst and second 

fnndamental forms of a hypersurface. We must then have: 

g%Kr,-i,j = = TriK^.,). 

We then deduce from (4.1) and (4.5) that there exists a constant C > 0 depending on Ai, 
i = 1,2, 3,4, 5, 6, 7, 8, and such that: 

(det^fjv)"^ < (det^°)"^exp(C't), (4.6) 


Now if we take in (4.1); 


{detg^)AQ > 1 , 


we will have for t sufficiently small; 

(det 5 f°)A 9 > exp(C't) > 1. 

This means, from (4.6) that, there exists ti > 0, such that, for 0 < f < fi, we have: 

(det^Ar)"^ < Ag. 

Then, using (4.3) and (4.7), we conclude that if T > 0 is such that: 


(4.7) 


0 < T < fi, BiT < Ai] i 


1,2,3,4,5, 6,7,8 


then we still have the inequalities (4.1) for n = N. Hence the iterated sequence (14) is nniformly 
bounded on [0, T[. □ 


4.3 Local existence 

Proposition 4.2. The initial values problem for the evolution system (2.51)...(2.59) has a 
unique local solution. 

Proof. Let [0, T[, T > 0, be the interval obtained in proposition 4.1. We prove that the iterated 
sequence 14 = {gn, K^, i4, E^, Un, 0n, Pn, Cn) converges uniformly on every [0, (5[c [0, T[, 5 > 0, 
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towards a solution V = (g, K, F, E, u, 0, p, e) of the evolution system. 


We take the difference between two consecutive terms of W, we use the fact that they have 
the same initial data and since the sequences > (Ki) are uniformly bounded we have, with 

(7 > 0 a constant: 


\gn+iit) - gn{t)\ -h \Kn+i{t) - 1 
+ \En+l{t) — En{t) \ + \Un+l{t) ' 
+ |Pn+l(0 ~ Pn{t)\ + |en-|-l(^) ~ 
+ \Kn{s)-Kn-l{s)\ + \Fn{s)- 
+ |'?in('S) ~ Un-l{s)\ -|- |0n('S) — 

For the same reasons we have: 


Ut)i 

+ 

iEn+l(t) 


- Un(t)j 

+ 

0n+l(^) 

-e„(«)| 

en(t)l 

< 

Cfo%ni 

~ 9n- 

■E^_,(s)j 

+ 

\En{s) - 

En-l{s) 

0n-l('S) 

+ 

|Pn(s) - 

Pn-l(s)| 


i(s)| (4.8) 

-f |e„(s) - en-i{s)\]ds. 


1 ^9n-\-l 

<^9n 1 _|_ 
dt 1 ' 

1 dKn-\-i 

. ^1 

F 

1 dF„+i 

_ dFn 1 

1 dt 

1 

dt 

dt ' 

1 dt 

dt 1 

1 1 f^-^n+1 

' 1 dt 

_ dEn 

+ 

1 dUn+1 

_ diLfi 1 

F 

1 <^©71+1 

_ d0n 1 

dt 

1 dt 

dt ' 

1 dt 

dt 1 

1 1 dpn-\-l 

' 1 dt 

_ dpn 1 

dt 1 

+ 1 

dSn^l 

dt 

_ den 1 
dt 1 

< 

Cl\9n 

9n—l 1 

F\K^- 

Kn-l\ 

+ 

\En- 

En-l\ 

+ 

\En — 

En-l\ 


+ \Un — Un-l\ + |0n ~ ©n-ll + \pn ~ Pn-l\ + l^n — Cn-ll]. 

For n G N, we set: 

Otnif) = \gn+l{t) - gn{t)\ + |^n+lW “ ^n(OI + l^n+lW “ ^n(f)| 

+ \En+l(t) — En(t)\ + — Un{t)\ + |0„+i(t) — Qn(t)\ (4-10) 

F \pn+l(t^ Pn(^)| F |Cn-|-l(f) 

(4.8) and (4.10) give: 


an{t) <C an-i{s)ds 


By induction on n > 2, we obtain, from(4.11) : 


\0in{t)\ < \\a2\ 


iCt) 


n—2 


'{n-2)\ 


< I|tt2| 




n—2 


(n-2)!’ 


for 0 < f < 5 and 0 < 5 < T. 

But the series ^ converges. Hence we obtain from (4.12) that: 


lim sup an(t) = 0. 

n->-|-oo o<t<5 


( 4 . 11 ) 


( 4 . 12 ) 
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Given the dehnition (4.10) of we conclude that every sequence {gn), {Kn), {Fn), {En), {un), 
(0n) {pn) and (e„) converges uniformly on every interval [0,5], 0 < h < T and we denote the 
different limits by g, K, F, E, u, 0, p and e which are continuous functions of t. 

Now from the inequality (4.9), we conclude similarly that the sequences of derivatives (^), 

(^) (^)’ (^) converge uniformly on [0,5], 0 < 5 < T. In 
these conditions, the functions g, K, E, E, u, 0, p and e are of class on [0,T[. Hence 
V = (g, K, E, E, u, 0, /9, e) is a local solution of the evolution system (2.51)...(2.59). 

We now prove that the solution is unique. 

Consider two solutions Vi and V 2 of the same initial values problem. Dehne a(t) = \ Vi{t) — V 2 {t)\ 
with Q!(0) = 0. Since the functions g, K, E, E, u, 0, p, e, (det ^r)"^ and ^ are bounded on [0, 5], 
0 < 5 < T, there exists a constant C > 0 such that: 

a{t) < (7 /" a{s)ds. 


Jo 

By Gronwall Lemma, we obtain a{t) = 0 since q;( 0) = 0. So Vi = V 2 and the local solution is 
unique. □ 


5 Global existence theorem 

We have to prove that, the solution V = {g, K, E, E, u, 0, p, e) and the functions 4 and (det g)~^ 
are bounded on every bounded interval. First of all, we prove the following important result on 
the mean curvature H = g^^Kij. 

Proposition 5.1. Let the function H = g^^Kij be bounded on [0,T*[ where T* < -|-cxd. Then 
the functions g, K, E, E,u,Q, p, e, (detp)“^, = a/1 + PijU^u^ and 4 are bounded on [0,T*[. 

Proof. We will use the following Lemma; 

Lemma 5.1. The mean curvature H = g^^Kij satisfies the following relation: 

= R + Ang^^Tij Tij) — 127r(roo + Foo) ~ 3A. (5-1) 

Proof of the Lemma 

Since F[ = g'^^Kij, the relation (2.17) gives: = 277*-^, from there we have: 
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then obtain (5.1) from equation (2.52) by a direct calculation using = 0, which implies 

= Too- 


Proof of Proposition 5.1 


Boundedness of l^r] on [0,r*[. 


— — A -|- 4:7ig^^{Tij -|- Tij) + 47r(roo + Poo)- 


By using the Hamiltonian constraint (2.61), (5.1) gives: 

dH 
dt 

But we have Tqo = g^^Tij > 0; Too > 0; g'^^Tij > 0. We then deduce from (5.2); 

^ - A. 

dt - ^ 

Integrating (5.3) on [0,f] for 0 < t < T* gives: 

H{t) > H{0) -At+ [ KijK^Us. 

Jo 

which shows, since H is bounded on [0,T*[ that we have: 

„T* 

J KijK^^ds < +00. 

Now the integration of equation (2.51) on [0;t], 0 < f < T*, gives: 

\g{t)\<\g{0)\+2 [ \K{s)\ds. 

Jo 

(5.6) gives, using (3.6), the inequality; 

\\gm<c 

(5.7), gives, using (3.7); 

\\gm<c 

By Gronwall Lemma, this implies: 

||9(()|| < Ci||9(0)||exp(c f\Kt,K->)ids I ,0 < « < T*; 


119(0)11+ / l|/f(9)IU.,ll9(9)ll* 

-'0 


|b(0)||+ / {K,,K^r^Us)\\ds 


This shows, using (5.5) that || 5 r||, and therefore \g\, are bounded on [0,T*[. 
• Boundedness of (det ^r)"^ on [0; T*[ 


The relation 


(5.2) 


(5.3) 


(5.4) 


(5.5) 


(5.6) 


(5.7) 


d 


— 


-[ln(det 9)]=9 
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gives, using equation (2.52) and H = 

^[In(det^)] = -2H. 

But \H\ is bounded on [0,T*[. Integrating on [0,f], 0 < t < T*[, there exists a constant C > 0 
such that; 

-C < ln(det g) <C 

from where we have: 

< det g < . 

Then det g and (det g)~^ are bounded on [0, T*[. 

• Boundedness of \K\ on [0,T*[ 

Since l^r] and (det g)~^ are bounded on [0, T*[, the expression (2.60) of Rij shows that R = g^^Rij 
is bounded on [0, T*[. We deduce from the Hamiltonian constraint (2.61), since Too > 0, Too > 0, 
that; 

R + H"^ -2A> KijK^A 

R and being bounded, so is KijK^K We then deduce from the inequality: 

\\K{s)\\<{K,,K^^)- 2 \\g{s)\\, 

that Ill'll and hence \K\ is bounded on [0,T*[, since and || 5 f|| are bounded. 

• Boundedness of T on [0,T*[ 

We use the Hamiltonian constraint (2.61) and Tqo > 0, Tqq > 0, to have: 

0 < max(167rroo; 167rToo) < 167r(roo + Too) = R + — 2A. 

This shows that Too and Too are bounded on [0,T*[. We deduce from (2.11) that; 

thus 

0 < ^gijE^E^ < Too, 

since Eij > 0. We use (3.9) to conclude that E is bounded, since {grsF ^^and \g\ are 
bounded. 
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• Boundedness of F on [0,T* 

Equation (2.53) gives, since l^r] and \E\ are bounded: 

l4'l [0,T*[. 

where (7 is a constant. Integrating on [0,t], 0 < t < T* gives the result. 

• Boundedness of p on [0,T*[ 

Since \E\, \g\, \K\, are bounded, equation (2.58) in p shows that, there exists two constants 

^4 > 0, I? > 0 such that; 

pit) < A [ p{s)ds + B, VtG[0,T*[. 

Jo 

By Gronwall Lemma, we have: 

pit) < < Be^^\ Vf e [0,T*[ 


and p is bounded. 


Boundedness of ^ and u* on [0,T*[ 


From (2.28) and (2.29), we have; 


p u 


A 1 


w 


p VP 4 


Since H = and ^ are bounded on [0; T*[, there exists a constant G > 0 such that; 




for 0 < f < r*. From where we deduce by integration; 


e ^ < pu° < on [0,T*[, 


which proves that; pu^ and ^ are bounded on [0,T*[. Now write, using (3.9): 
= from where we have: 


p|u*| < c(pM°)|p|2 
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and this prove that pu^ is bounded. Now multiply equation (2.55) in u* by p, and use the fact 
that pu^ is bounded to conclude that pii^ is also bounded. The relation 

u° = a/ 1 -h gijU^u^ 


gives, derivating and using equation (2.51); 

U* : 1 


u'^ 


pu^ 


which shows, since K, pu^, g, pii^ are bounded, that there exists a constant A > 0 such 
that; 

II 

< A. 




Integrating on [0,f], 0 < f < T* we obtain; 


1 < M°(t) < u/0)exp(ylT*). 

Hence is bounded on [0, T*\. Now write; 1 = ^ x nP and |u*| = y x to conclude that 1 

L ’ L p pU^ I I p 

and u* are bounded on [0,T*[. 

• Boundedness of on [0,T*[ 

Since \H\, p, \K\, (det g)~^ and u* are bounded on [0,T*[, by equations (2.56) and (2.57) in 
0°“, there exists two constant Ci and C 2 ' such that; 


|0°^| <C'i^|0°"| + C'2. 


q ;=0 

Summing in a and integrating on [0,t], 0 < t < T*, we obtain; 

3 3 pt ^ 

5^|0°“| < 5^|0°"(O)| + / ^|0°"|ds; 

q :=0 q :=0 q :=0 

and Gronwall Lemma gives: 

5^|0°“|(f) < [5^|0°“(O)| TdCsT*] exp(4C'iT*) on[0,T*[. 
Q=0 Va=o / 

which shows that each 0°" is bounded on [0,T*[. 

• Boundedness of e on [0, r*[ 


From the constraint (2.64) we have; 


e = 




which is bounded on [0,T*[, since and E are. This completes the proof of 
Proposition 5.1. 


□ 
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Theoreme 5.1. Suppose that A > 0 and that the mean curvature satisfies H{0) <0. Then the 
initial values problem for the coupled system of Einstein-Maxwell-pseudo tensor of pressure has 
a unique global solution. 


Proof. From Proposition 5.1, we have to prove that the mean curvature H is bounded on every 
interval [0,T*[ such that T* < -|-cxd. Consider the traceless tensor; 




(5.8) 


By a direct calculation, we have: 


- -Hf 

('J tj 2 

Using the hamiltonian constraint (2.61) and (5.9), we obtain; 

-2h = + 16n-(roo + Too) - R 


(5.9) 


(5.10) 


We have > 0, Too > 0, Too > 0. R.T. JANTZEN in [3] and R.WALD in [8] prove that, 
in the spacetimes considered here, we always have R < 0. We then deduce from (5.10); 

2 


this means; 


- 2A > 0, 
3 


> 3A. 


From there we have: 

H < —\/^, or H > \/^ 

but H is continuous and we have R(0) < 0, then we must have: 


H < -V^. 


Now if we take (5.4) in which we have > 0, we conclude that; 

Hit) > R(0) - At. 

Hence, on every interval [0,T*[, where T* < -|-cxd, we have: 

R(0) - AT* < H{f) < 

which prove that H is bounded on each interval [0, T*[. This ends the proof of theorem 5.1. □ 
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We end the paragraph by the following result; 

Proposition 5.2. For A < 0, there exists no global solution to the Einstein-Maxwell-pseudo- 
tensor of pressure initial values problem. 


Proof. Let A < 0 be given and suppose that the system has a global solution on the whole 
interval [0,-1-cxd[. 

By (5.3) and (5.9), we have: 


dH 

dt 


> + -H^ - A. 


But aija^^ > 0, then : 


dH 1 « 

— > -H^ -A. 
dt 3 

Since A < 0, we have —A > 0 and from (5.11) we deduce the inequalities: 

^ 

dt 3 


dH 

dt 


> -A. 


(5.13) gives, by integration on [0,f], t > 0; 


H{t) > H{0)-At. 


(5,11) 


(5.12) 

(5.13) 


(5,14) 


(5.14) shows that H{t) —> -|-cxd when t —> -|-cxd. So there exists to > 0 such that: H{to) > 0. 
Now we have, given (5.12); H{t) > W{t), for t > to, where W is every function satisfying; 



1 dt 3 

(5,15) 


\w{to) = H{to)>^. 

(5,16) 

A solution W of (5.15) on 

[fo,+oo[ is, given (5.16); 



W{f) = 

3-Hito)it-to) 

(5,17) 

(5.17) shows that: W{t) —)■ 

-|-cx) when t —)■< t* = to + > to. Hence H{t) - 

— > -|-CX) 


when t —)■< t*. This is impossible since the continuous function H is bounded on the compact 
set and cannot tend to -|-cxd when t tends to t*. So there can exist no global solution if 

A < 0. □ 
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6 Asymptotic behavior and geodesic completeness 

Here we extend the results of [4] to the charged case. We also study the asymptotic behavior 
of the charge e, the matter density p,the curvature tensor, the electromagnetic held and the 
pseudo-tensor of pressure. 

6.1 Asymptotic behavior 

Proposition 6.1. We have at late times: 


H{t) = + C>(e-2^‘) 

(6.1) 

H = 

(6,2) 

{detg)-^ = 0{e-^^^) 

(6,3) 

det g = 

(6,4) 

= 0(e-2T'‘) 

(6.5) 

a., = 

(6.6) 

roo(t) = 

(6,7) 

roo(f) = 0{e-^^^) 

(6,8) 

R = C>(e-2^‘) 

(6.9) 

KijK^^ = A + 

(6.10) 

g,,{t) = + 0(e-^‘)) 

(6.11) 

+ C>(e-^*)) 

(6.12) 

W^{t) = 

(6.13) 

EiE\t) = C>(e-2^*) 

(6.14) 

F,jE^^{t) = 0{e-^^^) 

(6.15) 

E\t) = 0(6^^*) 

(6.16) 
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F,, = G(e"^*) 

(6.17) 

is hounded 

(6.18) 

K,fit) = 0{e^^^) 

(6.19) 

p{u^f = C>(e-2^*) 

(6.20) 

p = 0(e-^'^*) 

(6.21) 

e = 0(e^^*) 

(6.22) 

e^j(t) = c>(e2(^-^)*) 

(6.23) 

nfit) = G(eS"*) 

(6.24) 

yfi- is bounded 

‘'J 

(6.25) 

000 = 000 = C>(e-2^*) 

(6.26) 

Tofit) = 0{e^^^) 

(6.27) 

To, = G(e®^*) 

(6.28) 

0Oj _ (7(e47t) 

(6.29) 

Tifit) = C>(e(®o'-2Dt) 

(6.30) 


where 7 ^ = 0 < // < 7 , {Gij) and are positive definite constant matrices. 


proof of ( 6 . 1 ) 

We have by (5.3) and (5.9), since aija^^ > 0; 

§>lW-A>0. (6.31) 


But; — A = — \/3A)(—H + \/^). We are in the case of global existence, then; 

H < —\/^. So; —H + ^/3K > 2^/^K. We then deduce from (6.31) that; 


We can write (6.32) as; 


where 


d(H + x/M) ^ + Tm) > 0; 

dt 


7 


2 


A 


(6.32) 

(6.33) 

(6.34) 
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Multiplying (6.33) by and integrate over [0,f]; f > 0, to obtain; 

> H{0) + 

which gives: 

\H + < |i7(0) + 

This proves (6.1) with 7 given by (6.34). □ 

Proof of (6.2) 

We have, for jj, g]0,7[: 

CLL (JjL (JjL 

which proves, given (6.1) and (6.35) since 0 < /r < 7 , that is integrable on 

conclude that —)■ 0 as f tends to -|-cxd and we obtain (6.2). □ 

Proof of (6.3) and (6.4) 


We have seen that; 


^[In (det^)] = -2H 


From H < —\/3A and (6.36) we have: 

67 < -2H < 67 + 2|if(0) + 


this means, given (6.37) that; 


67 < -^[In (det ^r)] < 67 - 1 - 2|i7(0) -|- v^^|e 

(JjL 


Integrating over [0,f], t > 0,we obtain: 


(det < det g < (det g^) exp ( — |i7(0) -|- \/^| ) 


7 


which proves (6.3) and (6.4). 


proof of (6.5) 

Since Tqo > 0, Tqo > 0, and —R > 0, (5.10) gives: 

-2A> > 0 . 

But we have, since H is increasing; 0 < —H + \/3A < —77(0) -|- y/SA. Then; 

0 < -772 -2A < 2| - 77(0) + V^\\H + V^\ 

3 


(6.35) 


(6.36) 
[0, -|-cxd[. We 

(6.37) 


(6.38) 


(6.39) 
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(6.5) follows then from (6.39) and (6.1). 
proof of (6.6) 

(6.6) is proven as in [4], using (6.5) and the notion of relative norms. 
proof of (6.7), (6.8) and (6.9) 

(5.10) gives, using (6.5); 

167r(roo + Too) - R + - 2A = ^(e-^^*) (6.40) 

But, Too > 0, Too > 0, < 0 and > 0. Hence (6.7), (6.8) and (6.9) follow from (6.40). 

proof of (6.10) 

(6.10) is a direct consequence of (5.9), (6.5) and (6.1). 
proof of (6.11) and (6.12) 

The proof of (6.11) and (6.12) are given by [4] using (6.6). 
proof of (6.13) 

We have; cr*-^ = (6.13) then follows from (6.12) and (6.6). 

proof of (6.14) and (6.15) 

We have; 

Too = g^^r,, = + \g^,E^E^ > 0 

Hence, 

0 < < Too , 0 < ]^gijE"E^ < Too- 

(6.14) and (6.15) are then consequences of (6.7). 
proof of (6.16) 

The second relation (3.9) gives, using (6.11) and (6.14); 

|T*| < {grsE^E^)^ \g\^ < [E^E^)^ \g\i < 

which proves (6.16). 
proof of (6.17) 

The integration of equation (2.53) gives; 

R,{t) = R,{0) + ^ gimE^{s)ds 

(6.17) then follows from (6.11) and (6.16) 
proof of (6.18) 
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Use (6-12) and (6.17). 

proof of (6.19) 

The relation (5.8) gives: 

I Fij I ^ I Cij I 2 I ^ 9ij 

(6.19) then follows from ( 6 . 6 ), (6.1) and (6.11). 
proof of ( 6 . 20 ) and ( 6 . 21 ) 

The relation (5.1) gives, since Too = g^^Tij: 




^2t,t 


dH 

dt 


-\- (3A — — R 127rToo -l- Stttoo 


The relations (6.1), (6.2), (6.7), ( 6 . 8 ), (6.9) and 0 < // < 7 , give: 


g^^R, = 0(e-2^‘). 


But we know that 


rr. 4 1 

Tij T '^pgiji 


thus: 

= ^g'^^UiUj + p> p{g"^UiUj -h 1) = p(m°)^ > P > 0. 

(6.41) gives (6.20) and (6.21). 
proof of ( 6 . 22 ) 

The constraint (2.64) i.e: Clj^E^ + evP = 0, (6.16) and > 1 give (6.22). 
proof of (6.23) 

(6.23) is given by (2.14) i.e: 0*^ = ^gij-, (6.11) and (6.21). 
proof of (6.24) 

From (2.48) we have: 


1 1 

Rj ~ — gikPjiE^E^ — -gijE^^Fik + gjiFi^F^^ 


this proves (6.24) using (6.11), (6.14), (6.16), (6.15), (6.17) and (6.18). 
proof of (6.25) 

It is a consequence of (2.19), using (6.11) and (6.12). 
proof of (6.26) 

From (2.12) we obtain: 

000 = Too — 2P('*^°)^- 


(6.41) 
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(6.26) is then a consequence of ( 6 . 8 ) and ( 6 . 20 ) since 0 < /x < 7 . 
proof of (6.27) 

From (2.48) we have; 

Toj = -E’^Fjk 

(6.27) is then a consequence of (6.16) and (6.17). 
proof of (6.28) 

We have, using the constraint (2.62); 


iTiToj = -STTToj + 


(6.28) is then a consequence of (6.27), (6.19), (6.24) and (6.12). 
proof of (6.29) 

From (6.12) we obtain; 


0 Oj _ _ -g^^Toi - 

(6.29) is then a consequence of (6.28), (6.12), (6.21) and (3.9). 
proof of (6.30) 

We have; 


nn 4 1 

Tij —pUiUj —ppij. 


From (3.9) we have; 

\Tij\ < cp{uy\gf + ^\pgij 


(6.30) then follows from (6.12), (6.20), (6.21), and (6.42) 
This ends the proof of proposition 6.1. 


(6.42) 


6.2 Geodesic completeness 


The geodesic equations for the metric (2.1) imply that along geodesics, the variables t, -u* 
satisfy a differential system which contains, between others, the equation; 

dt 


ds 


= u 


(6.43) 


where s is an affine parameter. The space-time will be geodesically complete if the affine pa¬ 
rameter s tends to -|-cxd as the time t tends to -|-cxd. Since (u°)^ = 1 -|- gijU^uf it will be enough 
if we prove that; 

rJM 1 

(6.44) 


— = (1 giju"u>) 2 > c> 0 . 
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since by integration we will have: s > cf -|- So, which proves that s tends to -|-cxd as t tends to 

-|-CX). 

We begin by the following result: 


Proposition 6.2. The quantity g^^UiUj satisfies the equation: 

[g^'^UiUj) = ^Hg^^UiUj + 


d / X 2 3eo 


Proof. From the evolution equation (2.51) and since gij is symmetric, we have: 

■^{g'^UiUj) = ^{gijU^'u^) = -2Kiju"u^ + 2gijiTuP 


3 /'* 1 

exp I - / r., . ds 1 E%i. 


(6.45) 


We can write this equation as: 


^{g^^UiUj) = -2K"^UiUj + 2iTui. 


(6.46) 


We have, using equation (2.55) of uh 


Ujj 


2Eu,, = 4K^^u,u. - 


u'^ 


2 pu' 


0 Wfe" 


2 p{u' 


0\2 jk-‘ 


(6.47) 


The last term is zero since EmivEv!' = 0, and the constraint (2.64) gives: 
Equation (6.47) then writes: 

3 e • 27b 

2Eui = AK^^UiU. - E^Ui --! . 

2 p 

We have from (6.46), (6.47) and (6.48): 

Jg^^u^ufi = 2K^^u,u, - ~E^u, - 


= —eu 


0 


(6.48) 


d , ■ 3 e • 2'y\,Eu^u, 

— (g^^Uiufi = 2K^^UiU. - E^Ui -- 1 

dV^ ^ 2p " uO 


(6.49) 


i . 


But we can write from expression (2.19) and 7^^ 




+gk„,{Cf;u^u^)u^ = 0 


since Off = —C^. Then, use the traceless tensor = Kij — -jgij, (6.49) gives 

— (g^^Uiufi = 2a^^UiUj H— Hg^^UiUj - E^Ui. 

dV ^ 3 ^ 2p 


(6.50) 
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Now we know by (6.26) that p satishes: 


P ( 4..n 




Integrating over [0,f], f > 0, we obtain; 


p{t) = p(0) exp - 


i /4 v 7 4„,a 


V„ 3 1 




+ 4 V ' * 


By (2.39) we have: 


Hence: 


e{t) = e(0) exp ( — 


4 4^0 




£W-££exp(5 T— 


(6.45) is then a direct consequence of (6.50) and (6.51). 

Proposition 6.3. If > | , the space-time is geodesically complete. 


(6.51) 

□ 


Proof. Since by (6.13) cr*-^ = 0{e ^^*), is bounded. The matrix in proposition 6.1 

being constant and positive dehnite, we have: 


< CG^^UiUj 


(6.52) 


where G is a constant. By (6.12) we can write: 


G^^u^Uj < Ge^^^g^^UiUj. 


(6.53) 


Since 5 ^ is a scalar product,we have: 

- E^Ui < {E^Ei) ^ {g^^UiUj) ^ . (6.54) 

So, (6.1), (6.45), (6.52), (6.53), (6.54) and (6.14) imply; 

^ {g'^UiUj) < (-27 + ce~‘^'^^)g^^UiUj + ce~'^^g"^UiUj + cexp (^-'jt + ^ ^ " • 

(6.55) 

From where, we deduce: 

^ [g^^UiUj) < (-27 + ce~'^^)g"^UiUj + cexp ^- 7 t j {g'^UiUj) ^ . (6.56) 

Sine > 1, we have ^ < t and - 7 t -h f ^ < t (-7 + f) 

So if; —7 -|- I < 0 or, if we set; cc = 7 — | > 0, we deduce from (6.56); 

— [g^^UiUj) < {—2u -f ce~^^)g^^UiUj -f [g^^UiUj) ^ (6.57) 
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Let us set: 

We have 

Hence, from (6.57) we obtain; 

^<ujZ+ (-2u + ce-^^)Z + 

at 

from where we have: 

dZ uj I uj I 1 . . 

— < ce 2 *Z-Fce 2 *Z 2. (6.58) 

dt 

But 

CO j. 1 UJ j. / UJ j. \ 1 r UJ j. u^ j. ~\ 

e-2*Z2 = e-4* (e-2*Z) < - [e"2* + e"2*Z] . 

We then deduced, from (6.58) that; 

+ <ce-^\Z + l). (6.59) 

dt 

(6.59) proves that Z = {g'^^UiUj) is bounded, and then, g'^^UiUj = gijU^u^ is bounded. So there 

exists H > 1 such that 

1 < < H. 

We then have (6.44) and this ends the proof of proposition 6.3. □ 

7 The positivity conditions 

Let ( "^X") be a future pointing vector. The quantity 

represents physically, the density of energy of a charged particle, measured by an observant 
whose speed is ( ^X"). Hence, this quantity should always be positive. Recall that a physical 
theory always has to satisfy at least one positivity condition [2]. There are three types of 
positivity conditions, they are, for ( ^X"), ( ^y“) future pointing vectors: 
a) The weak positivity condition which means: 

( + X/3) "x^ > 0 (7.1) 
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b) 


The strong positivity condition which means: 

‘‘.Y'’ > 0; 


(7,2) 


c) The dominant energy condition which means: 


( Vos + ‘y" > 0. 


(7,3) 


Obviously the dominant energy condition implies the weak energy condition, just setting 
4yo _ yy-g ]3y proving: 

Proposition 7.1. Let ( and ( be two future pointing vectors. Then: 

< 0. (7.4) 


Proof. Since ‘^Y^, the dehnition (2.1) of the metric ‘^g implies that (7.4) 

is equivalent to: 

- + gij ^Y^ < 0. (7.5) 

But ( ^X") and ( ^F") are future pointing, thus we have: 


< 0, > 0, ^F" ^F„ < 0, ^F° > 0 


or equivalently: 

0 < gij < ( ^X°)2; > 0; 0 < g^j ^F* ^F^' < ( ^F^)^; ^F° > 0. 

So, taking the square roots and the products: 

0 < {gij '^X^)^{gij ^F* ^F^')^ < ^F°. (7.6) 

Since g = {gij) is a scalar product, we have the Schwartz inequality: 

\g,, "X* *X^ < {g,, "X* ^X^)"^{g,, ^F* "F^)F (7.7) 

Hence, (7.4) is a direct consequence of (7.6) and (7.7). □ 

We now prove an important result. 

Proposition 7.2. Let ( "^X") and ( ^F“) be two future pointing vectors. The Maxwell tensor 
^Tag satisfies: 

^X“ ^F^ > 0. (7.8) 
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Proof. It shall be enough if we prove the inequality (7.8) in a particular frame. Guided by the 
electromagnetic held itself, we choose a frame of four future pointing vectors I = (/"); n = (n“); 
X = y = {y°‘) such that: 

= UafP = laX°‘ = la_y°‘ = = naX°‘ = 0. (7.9) 

The antisymmetric 2-form can be written in one of two forms; 

A B 

Fct/3 T "^{APyP ^yVa) ('^•^0) 

or 

C 

= -^{laXy - lyXa), (7.11) 

Where A, B, C are constants. It is important to choose the constants A, B, C such that; 

= -1; Xax'^ = = 1; x^y'^ = 0. (7.12) 

Let us consider the Maxwell tensor; 

+ "F„a "F/. (7.13) 

(7.11) gives, using (7.9) and (7.2); 

4 4 rp _ n-^TT' ^t-iA _ // 

^ ^ XfL U, -T Q,x r ^ ^ ^0:^/3 • 

(7.13) then gives: 

= ^IJp. (7.14) 

Let ( and ( ^F") be two future pointing vectors. We have: 

= j(L "X“)(/^ 

But since (/“), ( ^X“), ( ^Y°‘) are future pointing vectors, by (7.4) we have la < 0 and 
Ip < 0. Then ‘^Tap ^X°‘ ‘^Y^ > 0 and this proves proposition 7.2 □ 

Theoreme 7.1. The global solution of the Einstein-Maxwell system with pseudo-tensor of pres¬ 
sure satisfies: 

1°) The weak positivity condition; 

2°) The strong positivity condition if A> 
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Proof. - 

l'^) Let ( be a future pointing vector. By (7.8) we have: 

> 0. (7.15) 

Now by dehnition (2.12) of we have: 

‘^Ua ^U/3 + (7.16) 

If ( ^X") is a future pointing vector and since ( ^u") is by dehnition a future pointing vector, 
(7.4) implies, since p > 0: 

Qp s) ( V ( s "X^) > 0 (7.17) 

Now we use the dehnition of the pseudo-tensor of pressure ^©op. We have, using (2.13) and 
(2.14), for a fnture pointing vector: 

= ^©00 ( ^X^Y + 2 ^©oi ^X° ^X* + ^Qij ^X^ 

= ( ^ZoY ( ^X^Y + 2 (\ "Zo "X° "X* + ^pg,, "X* ^X^ 

= ( ^Z,Y ( ^X^Y + ( '^0 ^X* + ipp,, ^X* ^X^ 

= ( ^ZoY ( ^X^Y + ( ^^0 ^X“ - ^X°) + ^pp,, ^X' ^X^ 

o 

^©„p ^X“ *X^ = ( *Zo ^X°) ( ^X") + -ppij ^X^' (7.18) 

3 

But since ( "^Z") is future pointing, we have > 0, then "^Zq = — < 0, thus, *Zq "‘X^ < 

0 since ^X° > 0. 

Using once more (7.4) we have "^Zq. ^X" < 0, and hnally: ( "^Zq ^X°) ( "^Za ^X") > 0. Now the 
metric {gij) is positive dehnite, thus ^X® ^X-^ > 0. In conclusion (7.18) implies 
^©Q,p "‘X" *X^ > 0, and using (7.17) and (7.18) we obtain: 

Xp > 0. (7.19) 

Finally we have by (7.15) and (7.19): 

*x^ + *X^ = ( V„p + > 0 
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and the weak positivity condition is proved. 

2°) Let ( be a future pointing vector. 

The Einstein’s equations (2.8) imply: 

4^/5 + Stt ( 

By the weak positivity condition we know that; 

Stt ( ^X^ > 0. 


Let us see at what condition we will also have: 

4 ^^ > q 

By dehnition, we have: 

^R^p = ^g^^ "i?oo + "R^J = - ^Roo + 9^^ "R^r 


We know that ^Rij and Rij are linked by: 

^Ri, = R,, - dtK,j + HKij - 2KuK/. 
Contracting with g^^ we obtain; 


g^^ ^R,, = R- g^WtK,, + - 2W,,W*L 


Now we have, using equation (2.51); 


- K^^dtg^^ = dtH - 


From (7.24) and (7.25) we obtain; 

= R-dtH + H\ 


Now in the Einstein’s equations (2.8) take a = f3 = 0 to obtain; 

^Roo = — 2 —I" ^ + S'^( ^'^00 + ^^oo)- 
Then using (7.22), (7.26) and (7.27) we obtain: 

4/2 = !^_A_87r( Voo+ ^Too) + R-dtH + H\ 


(7.20) 


(7.21) 


(7.22) 


(7.23) 


(7.24) 


(7.25) 


(7.26) 


(7.27) 
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From where we deduce, since ‘^Tqo + Tqo > 0, R < 0 and dtH > 0 (see (6.31)); 


4 p 

Jl <-A+ H^ 


But H is increasing and negative, then < Ff^(O); hence: 

^ - A < -2A + (//(O))^ 

So we will have ^ — A < 0 if —2A-|- (i/(0))^ < 0 or A > Since ( "^X") is future pointing. 


we have by (7.4) ^Xq, < 0. In conclusion if A > 


MO)) 


■, then; 


— - A ] ^X“ ^X„ > 0 


(7.28) 


(7.28) and (7.21) imply 


^X" ^X^ > 0. 


This ends ends the proof of theorem (7.1)) 


□ 


8 Well-posedness 

According to Hadamard, a mathematical problem is well-posed if its solution exists, if the 
solution is unique, and if the solution is a continuous function of the initial data. 

We have only to prove the last point, in this paragraph we suppose eo > 0. The initial data is 
denoted Vq where: 

l/g = e D := 53(R) x 53(M) x ^ 3 (R) x R^°x]0, oo[x]0, oo[ 

where iS 3 (R) and Al 3 (R) are respectively the sets of 3 x 3 symmetric and antisymmetric matrices. 
We suppose that the initial data satisfy the constraints (2.61), (2.62), (2.63), (2.64). The global 
solution; 

A = (p,X,F,X,u\0O“,p,e) 

of the evolution system is also in D and is a function; 

S : D X [0, cxd[ —)■ D 

(Do,f) ^ SiV.R), 

We can also denote: 

SiVo,.) :[0,oo[ ^ D 

t ^ S{Vo,t), 

N.Noutchegueme and E.M.Zangue 36 





Dynamics of a plasma with pseudo-tensor of pressure 


such that; S'(Vo)O) = Vq. Every component of S is of class C^. Note that D is an open set of 
D = iS 3 (M) X iS 3 (M) X ^ 3 (M) X We can write: 

S :Q ^ CH[0,oo[,D) 

Do ^ 5(^0,.), 

If r is a compact set of [0, cx)[, we dehne the seminorm Pp on C^([0, cxd[, D) by; 

Pr{(p) = sup (I 99 I 1 ), (p E C^([0, cx)[, D) 
ter 

n ^ 

Where |IE|^ = ^ \Wi\ with W = {Wi,Wn) E C^([0, oo[, D) is a locally convex topological 

i=l 

vector space, whose topology r is generated by the family of seminorms Pp- We will prove that 
S is continuous from D to C^([0, cxd[, D), endowed with the topology r. Let p G N*. Set: 

Qp = B{0,2P)x]2-P,2P[x]2-P,2P[ 

where i3(0, 2^) is the ball of radius 2^ of iS 3 (M) x iS 3 (M) x ^ 3 (M) x Notice that: 

flp C Dp+i and U Dp — D 

pgn 

The function S will be continuous on D, if its restriction to every Dp is continuous. 

In what follows, we associate to every initial data the iterated sequence dehned in paragraph 
4. Recall that the iterated sequence converges to a unique solution if A > 0 and H{0) < 0. 

For every n eN, the iterated sequence 

Ln {,9ni P-m Fni Frit 

is a function; 

Vn : D X [0, cx)[ —)■ D 

{Vo,t) ^ Vn{Vo,t), 

Proposition 8.1. I 4 is continuons on Q x [0,cxd[. 

Proof. For Vo E fl, t E [0,cxd[, we have by dehnition Vo(ho,t) = Vq; hence Vq is continuous on 
D X [0, cxd[. Suppose that for n G N, Vq, Vi, •••, 14 are continuous, then the function; 

I4+i(Ro,t)= f foVM,s)ds + Vo (8.1) 

Jo 

where / dehnes the right hand side of the evolution system, is also continuous on D x [0, cx)[. 
Hence I 4 is continuous for every n G N. □ 
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Proposition 8.2. For p E N, 0 < T < -f cxd, the sequence (Vn)n is uniformly hounded on 
Vtp X [0, T[ 

Proof. Let us recall that; 

H{0)=g^’^^K^^ = Ho{Vo). ( 8 . 2 ) 

(8.2) shows that Hq is continuous on Qp x [0, T], hence it is bounded on it. Under the hypothesis 
that the iterated sequence satisfy the constraints (2.61), (2.62), (2.63), (2.64), we must have; 

Ho{Vo) < Hn{Vo,t) < E N, (UqU) e D x [0,oo[. 

Hq is bounded on Qp x [0,T[ and it does not depend on t. So there exists a constant Cp > 0 
depending only on p such that; 

—Cp < iLn(V(), f) < —\/ 3A, Vn G N, (Vo,t) E Dp x [OjU). 

This shoes that the sequence {Hn)nm is uniformly bounded on Dp x [o.ri. 

1°) boundedness of (5'n)nGN on Dp x [0,T[ 

The sequence (14)n satisfy the Hamiltonian constraint (2.61), that is; 

Hn — 167r(T„^00 + ^n,Oo) + — Rn- (8-3) 

From the fact that; A > 0, r„^oo ^ 0, T„,oo ^ 0, > 0, —Rn > 0 and {Hn)n 

bounded on Dp x [0,T[, (8.3) shows that the sequence {Kn^ijKfi)n of positive terms is 
bounded on Dp x [0,T[. Using the dehnition of the iterated sequence, the integration of 
equation (2.61) on [0,f], 0 < f < T gives; 

\9n+i{Vo,t)\<\9^\+2 [ \Kn{Vo,s)\ds. (8.4) 

Jo 

Using the notion of relative norm introduced in paragraph 3.3, we deduce from (8.4) the 
inequality; 

\gn+iiyo,t)\<C \g^\+ I \\gn{Vo,s)\\.\\Kn{Vo,s)\\g„(^Vo,s)ds (8.5) 

Jo 

where (7 > 0 is a constant. From (3.6) and (3.7), (8.5) implies; 

\\9n+i{Vo,t)\\<c\\\g^\\+ [\Kn,jKi^)^\9n{Vo,s)\\ds] . ( 8 . 6 ) 

Jo 
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But {Kn^ijKl^)n is uniformly bounded on Dp x [0,T[. Then (8.6) implies that there exists 
a constant Cp > 0, depending only on p, such that; 


\\gn+i{Vo,t)\ \ < Cp 


Ill7°ll+ / \\9n{Vo,s)\\ds 


, V(Do,t)eDpX[0,T[. 


By induction on n, we obtain; 


ll9„+i(u,«)ii<c,iiyi 

Finally, we obtain; 


1 -|- Cj)t -|- 


{Cpt)^ , , {Cptr , M„o,,(w 


+ ... + 


n 




n -I- 1 


V(Do,f) eDp 


\gn+i{Vo,t)\ < Cp\g^\{l + \g^\)exp{Cpt), Vn e N, {Vo,t) e Dp x [0,T[. (8.7) 


(8.7) shows that the sequence (g'n)n€N is uniformly bounded on Dp x 10, ri. 

2°) Boundedness of the sequence {Kn)n£f>i on Dp x [0,T[ 

Using once more the notion of relative norm, we have by (3.6); 

||^n(Uo,f)|| < (i^„,,,ir;-^')^||p„(UoT)||. (8.8) 

But {Kn,ijKl^)n and {gn)n are uniformly bounded on Dp x [0,T[. So is {Kn)n- 

3°) Boundedness of (detp^)ngN, ((det p^)"^)^ and {gl^)n on Dp x [0,r[ 

The relation (6.37) is equivalent to; 

det p = det exp(—2 f H{Vo,s)ds). 

Jo 

By dehnition of the iterated sequence, we must have; 

det5fn+i(Uo,f) = det5f°exp(-2 / Hn{Vo,s)ds), 

Jo 

From where, we deduce that, since {Hn)n is uniformly bounded on the compact Dp x [0, T], 
there exist a constant Cp^T > 0 depending only on p and T such that; 

det5(°exp(-C'p,T) < deign{Vo,t) < det exp(C'p,T), V {Vo,t) e Dp x [0,T[, 

This shows that (det gn)n and ((det gn)~^)n are uniformly bounded on Dp x [0, T[. we can 
conclude that (p)f )n is uniformly bounded on Dp x [0, T[. 


x|o.r|. 
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4°) Boundedness of {En)neN on Dp x [0,T[ 
From (8.3) we have: 


167r(rn,oo + Tnflo) — Rn — Kn^ijK'^ — 2A < — 2A. 


Since Kn,ijK^J_ > 0, i?„ < 0. But {Hn)n is uniformly bounded on Dp x [0, T[. Then so are 
(Tn,oo)n and {Tnfio)n- We know that; 

R,00 = + ^gnAjKEi 

then; 

0 < -gn^ijE^Eij^ < Tn^ij, 

since E^En^ij > 0. {gn,ijE'!^El)n is then uniformly bounded on Dp x [0,T[. Applying (3.9) 
and the fact that {gn,ijEl^E^)n and {gn)n are uniformly bounded on DpX [0, T[, we conclude 
that so is {En)n- 


5°) Boundedness of {Fn)n£f^ on Dp x [0, T[ 

By the dehnition of the iterated sequence, equation (2.53) implies: 

= Fj. + f C^^gn,ki{Vo,s)El{Vo,s)ds. (8.9) 

Jo 

But {gn)n and {En)n are uniformly bounded on Dp x [0, T[; by (8.9), (F„)n is also nniformly 
bounded on Dp x [0, T[. 


6 °) Boundedness of (u°)n, {j-)n, {pjn and «)„ on Dp x [0,T[ 


We saw that (2.28) and (2.29) give: 




p u 

and by integration on [0,t], 0 < t < T, we have: 

1 


pw = Pqw{0) exp 




H + Cl^-Ms 




By the dehnition of the iterated sequence, we must have; 


Pn+i^n+i = Po^^°(0)exp 


.3 1 


Ui 




( 8 . 10 ) 
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i ^ 

By (3.9) we have ^ < C \gn\‘^ ■, {gn)n and {Hn)n are uniformly bounded on Dp x P.ri, 
and > 1. So, from (8.10) we deduce that there exist a constant Cp^T > 0 such that: 

poM°(0) exp(-C'p,r) < pn+iK+i < PoW°(0) exp(C'p,T), on Dp x [0, T[. 


So {pnuDn and 


' Pn'^ri 


are uniformly bounded on Dp x [0, T[. But > 1, then {pn)n is 


also uniformly bounded on Dp x [0,T[. Now we have by (3.9): 

PnU'- 


u 




PnU^ 


<C\gr, 


But {pnU^)n and {gn)n are uniformly bounded on Dp x [0, T[. Hence so is {pnu\^n- Now 
consider the equation in deduced from (2.55). We multiply this equation by pn and 
conclude that the sequence {pnuh+i)n is uniformly bounded on Dp x [0,T[. Now derive 
the relation = a/1 + gijU^u^ and obtain, using (2.51): 


iP u* pii^ 


From there, we deduce that for the iterated sequence, we must have: 


u. 


'n+l 


rilJ 


= -K 


u. 


n+l 


, (Pn<) < 


Since (iFn)n, (^)n, {gn)n, (PnOn, {-^)n are Uniformly bounded on Dp x [0,T[ there 
exists a constant Cp^T > 0 such that: 


u. 


n+l 


U. 


n+l 


— 


From there, we deduce that: 

1 < u^n+i — exp(C'p^'rT), on Dp x [0,T[ 

then {u^j^Pn is uniformly bounded on Dp x [0,T[. Now -y = X thus (^)n is 
uniformly bounded on Dp x [0,T[. We have ^ x and are nniformly 

'^n '^n 

bounded on Dp x [0,T[. Then (u)j)n is uniformly bounded on Dp x 10, T|. 

In conclusion, the sequences {pn)n, (^)n, and (u)j)n are uniformly bounded on 

Dp X [0, T[. 

7°) Boundedness of (e„)„gN on Dp x [0,T[ 

The sequence (e„)„ from the iterated sequence satishes constraint (2.64), hence: 

1 

„ _r<i j+k 

On — Q 

“n 

Since > 1 and (i7n)n is uniformly bounded on Dp x [0,T[, so is (e„)„. 
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8 °) Boundedness of on x [0)^[ 

By the equation dehning the iterated sequence, since {Hn)n, iPn)n, (n°)n, {gn)n, {Kn^ij)n, 
inn,ij)n and {u\)n are uniformly bounded on Dp x [ 0 , r[, we can conclude that, there exist 
two constants ^p,r > 0 , and -Bp,r > 0 such that; 


lCi(^o,f)| < |0°°(Do,O)|+Ap,r^*^|e(K),5)|ds + i?p,T onDpX [0,T[ 

and 

lCi(^o,i)l < l0°*(Do,O)l+Ap,T /V|C(K),5)lds + 5p,T onDpX [0,T[ 

>^0 „=o 

From where we deduce that: 

5hCi(U,«)| <4/lp,T f‘j2K(Vi>,s)\ds + 4B,,T + J2\e^(Vi>,0)\ wiOpX 10,T[ 

By induction on n, we conclude that; 

5^|Ci(H,f)|< (4ylp,T + 5^|0°“(Do,O)|+4Sp,T)exp(4TC'p,T), onDpX[0,T[. 

0=0 V 0=0 / 

(841) 

Where Cp^T is a constant. Now we know that 

q :=0 


and |Vo|i is continuous on the compact set Dp x [0,T]; then it is bounded on it. In 
conclusion, ( 8 . 11 ) shows that every sequence is uniformly bounded on Dp x |0,T|. 

This ends the proof of proposition 8.2 


□ 

Theoreme 8.1. The Cauchy problem for the Einstein-Maxwell system with pseudo-tensor of 
pressure is well-posed in the sense of Hadamard. 

Proof. We prove that the function; 

S :Q ^ Ci([0,oo[,D) 

Do ^ ^(Do, . ), 
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is continuous. We know that it will be enough if we prove that its restriction to every Dp is 
continuous. 

Since > 1 and since the sequences (14)^ ( —) are uniformly bounded on Dp x P.ri, 

/ n 

by the dehnition of the iterated sequence, there exist a constant Cp^x > 0 such that; 


We have: 


— 


t) - ^, V Do, Do e Dp, n e N, f e [0, T[. 


Dn(Do,0) = Do ; D4Do,0) = Do. 
Then, integrating this inequation on [0,f], 0 < f < T, we have; 

|D„+i(Do,t)-D„+l(DoT)|l<4p,T [ Vn{Vo,s)-Vn{Vo,s) 

Jo 

By induction on n, we deduce from this inequation that; 

{tCp^T 


ds + 


Do-Do 


|D„(Do, t) - D„(Do, t)|i < |Do - Doll 


i=0 


from where we deduce: 


|D„+i(Do,t) - D„+i(Do,f)|i < |Do - Do|iexp(TC'p,T) 

Taking the limit when n —)■ -|-cxd, we obtain; 

\S{Vo,t) - S{Vo,t)\, < |Do-Do|iexp(TC'p,T),V Do, Do G Dp ,f G [0,T], 

Using the seminorm P[o,r] we obtain; 

P[0;T]{S{Vo, .) - >5(Do, .)) < |Do - Do|iexp(TC'p,T),V Do, Do e Dp. (8.12) 

which proves that Vq i —> 5'(D, •) is continuous from Dp to U^([0, cxd[, D). 

This proves that the Cauchy problem for the Einstein-Maxwell system with pseudo-tensor of 
pressure is well-posed. □ 
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Conclusion and fnture investigation 


We have proved global existence of solutions to the evolution system [S) when the cosmological 
constant A is strictly positive and H{0) < 0. We have also proved the geodesic completeness 
(when > |), and determined the asymptotic behavior (of the spacetimes) in the neighbor¬ 
hood of -l-cxD in the case of global existence. On the other hand, we have hnally proved that the 
problem is well posed in the sense of Hadamard. 
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APPENDIX 
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Point 1 of Proposition 3.1 

1) Evolution of ^4 = i? — KijK^^ + — 167r(roo + Poo) — 2A 

1.1) Evolution of H 

From definition (2.18) of iP, we have: 

dtH = 

But equation (2.52) gives: 

= g^^ [R,, + HK,, - 2K]Ku - + R,) + 47r(-Too + g'^Tlm)g^j - 

= R + H^- 2K^^Ka - Sng^^r^j - ^ng^^R, + 127r(-Too + g^'^Ti^) - 3A. 

From the two relations, we have the evolution of H: 

dtH = R + + Airg^^Tij — Stttoo — 127rToo — 3A. (^-1) 

1.2) Evolution of Tqo + Poo 


We use the index ( ^) to write the conservation laws: 

"V„( + "P“^) = ^g^^ V" + ^Tx,) = 0. 

Now, using (2.7) we have: 

Then use (2.2), (2.3) and (2.19) to obtain the evolution of Tqo + Poo^ 

^(roo + Poo) = H{roo + Poo) - g'^7-j{rok + Pofc) + + P.,) (A.2). 

1.3) Evolution of R 

We have by dehnition: 


1 


R 

^R 


X 

a, 06 


g^^Rij 

- ^esi ^ 7 ^„) + ^ 7 ^„ - ^ 7 ^ (^- 3 ) 
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Where are dehned by; 


1 _ 4^A 


ex 


{AA) 


with; 


^0 _ 4^A _ 4/^A _ ri. 4/^fc _ 

'-'a/3 “ '-'03 “ '-'aO “ '-'ii “ 


''03 ~ '-'aO 




(AS) 


In particular, we deduce from (A.3); 


4 

J,i0 


“ei/S'o)- ‘e„(‘‘7',)+ S'PTJi- “tIS;-- ‘CJ-T. 

-|(7t) - (A-yL - 


TJ 


Hence; 




Now we have from the Codazzi relation; 


^Rxi,ji iKij 'VjKii. 


But here we have n = dp, this implies; 


Roi,jl — ~'^lKij + jKil 


(4.6) 


or; 

X',o = -V%' + V,W'; 

(A.6) and (A.7) give; 

1(7^ = v'A'„ - V,A''- V.A'' 

Now a direct calculation gives; 

Vz ( Ay; 


dt 


^Rij. 
dt ^ 


{A.7) 


(A8). 


We deduce from (A.8) that; 


d 


-R,, = ViVKij - ViVjKl - ViViK] 
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We deduce: 


and; 


But; 


= -2g"' [7,TV‘A',„„ - 7taV’"A7„ - 7LV’’‘A„.] 
= -2 [7,TVAt - 7 LV’”a 1 - g‘"7LV’"A„J 

= 2g»g^^V‘K,t. 


R = g‘>R. 


V 


dtR = 2K^^Rig + g^^—Rig 

ij y dt 


and the evolution of R: 


d,R = 2K»Rij + 29«7‘V'An. 


(yl.9) 


1.4) Evolution of 
We have; 


dt{Kg^K^^) 




+ Kg, 


K^^Kij + K, 


Hence: 

dtiKgjK^^) = 2K^^ki, + AKg^K^^Kl 

Then equation (2.52) gives the evolution of KgjK'^^: 


dt{Kg,K^^)=2K^mg,+2HKg,K^^-167rK^kRj+Tgj)+S7rH{-Too+g^^Tim)-2HA. 


1.5) Evolution of A 
We have: 


dtA = dtR - dt{Kg,K^^) + 2HdtH - Wndtiroo + Too) 


and using (A.9), (A.10), (A.l) and (A.2), we obtain the evolution of A: 


dtA = 2HA + 2g‘’j^gAt. 


(All) 


(AlO) 
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2) Evolution of Aj = V^Kij -|- 87r(roj -|- Tqj) 

We have to compute: 

^ , and {toj + Toj). 

If we take /d = j in the conservation laws (2.21), we obtain; 

Va( = ^Vo( + ^r^') = 0, 

or, using (2.7); 

+ ^T°^') + ^7o\( 

+ Hxi + ^d( + ^r^') + V' 

+ Hx( + ^r^) 

= 0 . 

From there we deduce, using (2.19) that: 

^(^Oi ^ ^ ^ ^ 2K/(r°* + T°*). 

Now; 

Ton + To„ = ^goa "^n/3( -h "T"^) = -^n,(r°^' + T°^) 

Then, using yfi'nj = —‘^Knj given by equation (2.51), we obtain; 

^(ron + Ton) = 2Knj{T^^ + T°^) - + T^n. 

We then use (A. 12) to obtain; 

+ Toj) = H{Toj + Toj) + V*(rjj -|- Tij). (A.13) 

Now we have by Bianchi identities: 

^ = 0 . 

But using (2.19), with /? = j, we have: 

4v/ 3 4^ _ 4y/3 4^a _ 4y^ 4^a _ 4^*j [ 4^a 4^A _ 4^A^ 4^aj _ , 


j _|_ irpXj^ 


(4.12) 
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Then we have: 


But; 


= 0 . 


= 0 ) ^ = 0 ), 

and developping, we obtain; 

- Ki = 0. (A.14) 

This means that we have to compute 

and 


We know that: 


4 


R, 


'a0 




^e/3 (SL) + ^7^ - ^7^^^7L - 


A 

lia' 


We then obtain, using once more (2.19); 


(A. 


4^0i _ 4^0a 4^i/3 4^^^ 

= R[KriL-Ri"iS 


Thus; 

4i?oi = (A. 16) 

If we consider formula (A.15), we deduce using (2.19) that; 

= -9“g^ 'a, A-„ + HK‘i + g‘V' [ S?, - Sf, - CJ, S«J - 2A-yff» 

i.e., using (2.60); 

4^0 _ ^ - g^^g^^dtKki. (A.17) 
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We then deduce, using (A. 16); 

= -H *R^ - Ki + ^7*;- ^R^ + ^R^ 

and we deduce from (A. 17) and (2.19); 

= -KjWkR^ + V,R^ - 2V,{K;KR - g^Pg^'^V^idtKpg). (A.18) 

We obtain from (A.14), (A.16) and (A.18); 

- KlWkR^ - KiVkR’^ + V^R^ - 2V^{KIKR - g^Pg^^V^{d,K„) = 0 
We now use equation (2.52) to express dtKpq and we obtain; 

j^ig^’^RKR - 2KlVkR^ + V^R^ - 2V^{KIKR - 

g^^g^^^V, [Rpq + HKpq - 2K^Kpk - 8n{Tpq + Rq) + 47r(-Too + g^^RRg^q - kg^q] = 0 

Or, since the functions depend only on t, and -^g^^ = 2K^^: 

g^^j^{RKR+2R^RRu-2KlVkR^ + -2V^{K;kR 

-ViR^ - + 2Vr{R^Ki) + 8TiVi{R + R) = 0 

or; 

9jn g^’^j^RKR - HVR^R87rV,{R + rn =0, 

that is; 

= HRRj - 8TiV\Tiq + R). (A.19) 

We hnally obtain, using (A. 13) and (A. 19); 

4.4, = ^(V/f,)4 84(r„,4T„,) 

= HW^Rq - 8TiW\nq + R) + 87r [7/(ro, + R) + W\nq + R)] 

= H [V*A"jj + 87r(roj + TQj)~\ ■ 

and the evolution of Aj: 

4-A, = HA. 

dt ^ ^ 
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3) Evolution of Ajfc = + C^j^Fu + Cl^Fji 

Using equation (2.53) in Fij, we have: 

= Cld,(F,,) + C>A{F,,) + Cld,(F,i) 

= c‘, (CZar^iE') + {C",'g,^,E‘) + CL (qUj,n,E‘) 

_ / yom I yom I \ jpl 

— “T ^jk^il “T ^ki^jk) yml^ 

= 0 , 


from the equality of Jacobi. Hence: 


dtAijk = 0. 

4) Evolution of H = Cl^F^^ + evP 


Set 


= ClkF^^ -1- eu“ 


We have = B. Setting = eu”, we have from equation (2.9): 


= CyVa^F 


= CIJK 


Thus: 

or: 

jB+ = c‘qE 

i.e. 


(It 


B + 


"ihB + H 




B^ 


ClkJ 


k. 


Use (2.20) to obtain: 


d 

-B-HB 

dt 


q[qp'‘ + eu=]=c‘,eu'‘. 


N.Noutchegueme and E.M.Zangue 


52 





Dynamics of a plasma with pseudo-tensor of pressure 


or; 


±B-HB + Cifil.Ff + CifM’ = CieM*. 

But 




Hence: 


and the evolution of B; 


This ends the appendix. □ 


-HB=0 
dt 



HB. 
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